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We study an opinion dynamics model in which agents reach compromise via pairwise interactions. 
When the opinions of two agents are sufficiently close, they both acquire the average of their initial 
opinions; otherwise, they do not interact. Generically, the system reaches a steady state with a 
finite number of isolated, noninteracting opinion clusters ("parties"). As the initial opinion range 
increases, the number of such parties undergoes a periodic sequence of bifurcations. Both major and 
minor parties emerge, and these are organized in alternating pattern. This behavior is illuminated 
by considering discrete opinion states. 
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I. INTRODUCTION 

In a society, people typically have a wide range of 
opinions. However, individual opinions on a particular 
issue are not static but rather evolve due to the influ- 
ences of acquaintances or other external factors. In prin- 
ciple, opinions could evolve ad infinitum, consensus could 
emerge, or a population could reach a state that consists 
of a finite set of distinct opinion clusters, or "parties". 

It is natural to discuss this process within the frame- 
work of interacting particle systems [Q, • A classic ex- 
ample is the voter model where agents, who possess two 
possible opinions, adopt the state of a randomly selected 
neighbor ]j| ||. Individual opinions evolve until con- 
sensus is eventually reached, and the probability that 
a given opinion ultimately wins is equal to the initial 
fraction of agents with that opinion |]. Several other 
Ising-type opinion models, incorporating more realistic 
features, have been proposed recently P, Ell m, |l 5 |l0| . 

In this paper, we study the compromise model, a sim- 
ple model for the evolution of opinions in a heterogeneous 
population jll], [l^] . To account for the diversity of the 
population, the opinion is either a real-valued variable 
or a discrete variable with many states |l3|, [l4|, [Tq ]. To 
mimic the natural human tendency for reaching a fair 
compromise, in an interaction between two agents, both 
acquire the average of their initial opinions. Last, to in- 
corporate self confidence or conviction in one's own opin- 
ion, interactions between agents whose opinion difference 
is larger than some threshold are forbidden. 

Monte Carlo simulations of the compromise model 
have shown that either consensus or diversity can arise, 
depending on system parameters [jll], Q. Here, we inves- 



tigate the compromise model using numerical integration 
of the governing rate equations for continuum opinions 
and analytical solutions for discrete opinions. Numerical 
integration is more efficient than direct numerical simula- 
tion and provides better resolution of the time dependent 
and steady state behaviors. 

We find that the compromise model exhibits a rich be- 
havior. In the long-time limit, the system condenses into 
a finite set of equally-spaced opinion clusters (parties), 
with the population in adjacent clusters alternating be- 
tween two values that differ by 4 orders of magnitude. 
As the initial range of opinions grows, the number of 
parties increases via a periodic sequence of bifurcations. 
The corresponding period governs the basic features of 
the emergent structure, namely, the size of the major 
clusters, and their separation. Near bifurcation points, 
the size of minor clusters vanishes algebraically, and we 
provide a heuristic explanation for this behavior. 

Underlying the compromise model is a stochastic aver- 
aging process. Closely related averaging processes natu- 
rally arise in diverse systems, including one-dimensional 
inelastic collisions [[l6[ [l7j , dynamics of headways in traf- 
fic flows m |, mass transport pofl , force fluctuations 

in bead packs ]2l[, wealth exchange processes |2^, p3[ , 
and the Hammersley process [^4[ |2f| . While our findings 
are discussed in the framework of opinion dynamics, they 
may very well be relevant in these different contexts. 
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In Sec. II, we describe the numerical integration of the 
rate equations for the opinion probability density and the 
resulting bifurcations. In Sec. Ill, we examine systems 
with a finite number N of discrete and equally-spaced 
opinions. When N is relatively small, these systems can 
be treated analytically, thereby illuminating the behavior 
in the continuum case. Generally, consensus is reached 
for small enough N, while a state with several distinct 
non-interacting clusters is reached for large N. We con- 
clude in Sec. IV. 
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II. THE CONTINUUM VERSION 

In the continuum version of the compromise model, 
each agent is initially assigned an opinion x from some 
specified distribution. Randomly selected pairs of agents 
undergo sequential interactions. Such interactions are 
restricted to agents whose opinion difference lies below 
a threshold that is set to unity without loss of general- 
ity. When agents with opinions X\ and X2 interact, both 
acquire the average opinion: 



{x\,x 2 ) 



xi + x 2 xi + x 2 



\x 2 -xx\<l, (1) 



while if \x 2 — x\\ > 1, no interaction occurs. This model 
is essentially identical to that of Refs. 0, [L2| . 

Let us denote by P(x, t) dx the fraction of agents that 
have opinions in the range [x, x + dx] at time t. The 
distribution P(x,t) evolves according to the rate equation 



J^-POM) = / / dxidx 2 P{x u t)P{x 2 , 



t) 



2:2 1 < 1 
xi + x 2 



5 x 



5(x — x\) 



(2) 



The quadratic integrand reflects the binary nature of the 
interaction and the gain and loss terms reflect the pro- 
cess (Q). This basic dynamical rule conserves the total 
mass and the mean opinion. That is, M$ and M\, the 
first two moments of the opinion distribution are con- 
served, where Mk (t) = J dxx k P{x, t) is the fcth moment 
of the distribution. We restrict our attention to flat ini- 
tial distributions Po(x) = P(x,0) = 1 for x E [-A, A]. 
Our goal is to determine the nature of the final state 
Poo(x) = P{x, 00). 

When all agents interact, namely, when A < 1/2, the 
rate equations are integrable JToj , |l7| . In particular, the 
second moment obeys M 2 + MqM 2 /2 = Mj 2 , where the 
overdot denotes time derivative. Using M\ = 0, we find 
that the second moment vanishes exponentially in time: 



M 2 (t) = Af 2 (0)e- Mot / 2 , 



(3) 



with Mq = 2A. Thus all agents approach the center 
opinion and the system eventually reaches the consensus 



Poo Or) =M S(x) 



(4) 



The distribution P(x,t) approaches the localized state 



in a self-similar fashion, P{x, t) 

T 1 / 2 , 



2M> 



1 



with variance w — M 2 /Mq and scaling variable z = 
x/w§Pj§. 

For larger values of A, the opinion distribution does 
not condense into a single cluster, but rather the dis- 
tribution evolves into "patches" that are separated by a 
distance larger than one. This behavior results from an 
instability that propagates from the boundary toward the 
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FIG. 1: Evolution of the opinion distribution for A = 4.3 
where four major clusters ultimately arise (see Fig. 2). Shown 
is P(x, t) versus x for times t = 0.5 (bottom), 3, 6, and 9 (top). 



center (Fig. |l|). Once each patch is isolated, it then sep- 
arately evolves into a delta function as in the A < 1/2 
case. The final distribution consists of a series of non- 
interacting clusters at locations Xi with masses mf. 



Poo(x) 



rrii 5{x — Xi) 



(5) 



with ^2 mi = Mq = 2A and m>iXi — Mi = to satisfy 
the conservation laws. 

Our goal is to understand basic characteristics of the 
final state. How many clusters arise? Where are they 
located (in opinion space)? What are their masses? As 
we shall see, the answers to these questions depend in 
a surprisingly complex manner on the single control pa- 
rameter, the initial opinion range A. 



A. Cluster Locations 

To determine how the final state depends on A, we nu- 
merically integrated the rate equation (|^) by discretizing 
x into 400A equally-spaced states. The range < A < 10 
was investigated using a fine mesh (0.0025 increments). 
The rate equations were integrated using a fourth or- 
der Adams-Bashforth method |^6| up to a sufficiently 
long time that the probability distribution separated into 
nonintcracting patches. Then, the two conservation laws 
were invoked to determine the ultimate mass and loca- 
tion of each patch. The accuracy was 10 -9 in P(x,t). 

The cluster locations exhibit a striking regularity, as 
seen in plotting Xi versus A (Fig. |^). There are three 
types of clusters: major clusters (mass M > 1), minor 
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FIG. 2: Location of final state clusters versus the initial opin- 
ion range A. The three types of clusters, defined in the text, 
are noted. The vertical arrows indicate the location of the 
first 3 bifurcations. 

clusters (mass m < 10~ 2 ), and a central cluster located 
exactly at x = 0. The number of clusters grows via a 
series of bifurcations. When A < 1/2, the final state 
is a single peak located at the origin, and this situation 
persists as long as A < 1. When A exceeds one, two new 
clusters are born at the extreme edges, x » ±A. As A 
increases, further bifurcations of three basic types occur: 

1. Nucleation of a symmetric pair of clusters: — > 
{—x, x} with i=l. 

2. Annihilation of a central cluster and simultaneous 
nucleation of a symmetric pair of clusters: {0} — > 
{—x, x} with x « 0.75. 

3. Nucleation of a central cluster: — * {0}. 

The bifurcations always occur in a periodic order: 
1,2, 3, 1,2, 3,.... Numerically, the first 4 generations of 
bifurcations are located at 



(Ai,A 2 ,A 3 



(1.000,1.871,2.248), 
(3.289,4.079,4.455), 
(5.496,6.259,6.638), 
(7.676,8.431,8.810). 



Successive bifurcation points of the same type are all sep- 
arated by the same distance: Ai(n + 1) — Aj(n) — > const. 
Also, the distance between different types of bifurcations 
within the same generation eventually becomes constant. 
Thus the bifurcation diagram, with all its intricate fea- 
tures, repeats in a periodic manner 



x(A) = x(A + L), 



(6) 



with period L rj 2.155. The period was estimated by 
extrapolating the differences in the locations of the first 
few transitions to oo. 
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FIG. 3: The masses of the final clusters versus their location 
for the cases of A w 7.8 and A + L ~ 10. The clusters in the 
two systems coincide, except that the larger system has four 
more clusters. The scale of the inter-cluster separations and 
the masses of the major clusters are indicated. 



In type-1 (2) bifurcations, branches of minor (major) 
clusters nucleate near the origin, and these persist for all 
larger A. As each branch evolves, notice that it exhibits 
a large curvature change or a kink due to the effect of a 
subsequent bifurcation. For |x| ^ 2, the branch growth 
is practically linear with a slope commensurate with the 
opinion range: \dx/dA\ — > 1. 

The periodic behavior further implies that the separa- 
tion between clusters becomes constant. Moreover, when 
a system of size A is compared with a system of size 
A + L, cluster locations in the smaller system coincide 
with the larger one, as shown in Fig. ||. The larger sys- 
tem, however, contains two additional pairs of major and 
minor clusters. Thus, the period L governs the overall 
number of clusters and the separation between them. For 
A ^> 1 there are 4A/L clusters, with neighboring clusters 
separated (approximately) by distance L/2. 



B. Cluster Masses 

An even richer picture emerges when the cluster masses 
are considered. First, the cluster masses vary periodically 
in the initial opinion range, that is, m(A) = m(A + L), 
as seen in Fig. |4|. Second, clusters are organized in an 
alternating major-minor pattern (see Figs, g an d ||). For 
large A, each cluster mass approaches a constant value. 
The major clusters, which contain nearly the entire mass 
in the system, saturate at a value equal to the period, 
M — » L. The masses of the minor clusters approach a 
much smaller level: m — ► 3 • 10~ 4 (see Fig. [|). This 
minute mass implies that a sufficiently large population 
is needed for minor clusters to exist. 

The central cluster is special. Its mass never becomes 
constant but instead varies in a periodic manner with A 
(Fig. ||) . A central cluster nucleates with an infinitesimal 
mass at a type-3 bifurcation, grows slowly for a while, 
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FIG. 4: Cluster mass versus opinion range. The central clus- 
ters (periodic variation) and the major cluster are shown on a 
linear scale (top), while minor and central clusters are shown 
on a logarithmic scale (bottom). 



then it undergoes an explosive growth until its mass be- 
comes of order unity. Finally, its mass grows linearly 
with A. At some threshold, the central cluster splits into 
two major clusters via a type-2 bifurcation (Fig. [|). This 
birth-and-death pattern repeats ad infinitum. 

The minor clusters exhibit two subtle features. First, 
the mass of the most extreme cluster saturates to a mass 
m! that is approximately one order of magnitude greater 
than all other minor clusters. Second, the mass of the 
minor clusters varies non-monotonically with A and there 
is a small range of A, where the mass of a newly-born 
minor cluster suddenly drops (Fig. |j) before the mass 
saturates to a constant value. We are unable to resolve 
whether there is a finite gap or just a singular point where 
the mass vanishes. 

At type-1 and type-3 bifurcations, new clusters form, 
and the mass of these nascent clusters vanishes alge- 
braically according to 



(A - A„) £ 



(7) 



as A — > A„. The exponent depends only on the type n 
of the bifurcation point; numerically we find a,\ « 3 and 
c*3 w 4 (Fig. ||). We now give a heuristic explanation for 
this behavior. 

To understand the behavior near a type-1 bifurcation, 
consider the very first one at Ai = 1. Let A = 1 + e 
with e — > 0. It is convenient to divide the total opinion 
range (—A, A) into a central subinterval (—1, 1) and two 
boundary subintervals: (1, 1 + e) and (—1 — e, —1). Let 
m(t) be the mass in a boundary subinterval. Initially, 
m(0) = e. Such mass is lost due to interaction with one 
half of the central subinterval. As a result, ra = — m, 
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FIG. 5: Critical behavior for the masses of the minor clus- 
ters at type-1 (top) and type-3 (bottom) bifurcations. The 
straight lines have slopes 3 and 4. 



which gives 



m{t) 



e e 



(8) 



On the other hand, the mass of the central subinter- 
val gets concentrated in a region near the origin whose 
spread w(t) decreases with time. At some moment tf 
the separation between masses in the central and bound- 
ary subintervals exceeds unity. We anticipate that the 
mass in the boundary subinterval converges to its center 
x = 1 + e/2, and hence, this critical separation occurs 
when w(tf) ~ e/2. For t 3> tf, the interaction between 
the two subintervals stops and the mass of the emerging 
minor cluster freezes at m/ ~ ee _t/ . 

The spread w(t) can be estimated by noting that to 
zcroth order in e, (i) the central subinterval is not affected 
by the boundary subintervals, and (ii) eventually, almost 
all agents are within the interaction range. Therefore, the 
asymptotic behavior is the same as in the case A < 1/2, 
and the spread follows directly from the second moment 



w(t) ~ M. 



1/2 



"*/ 2 since A 



the stopping criteria, w(tf) 



1 + e = 1. Using 
/ 2 ^ £j the final minor 
cluster mass is m/ ~ ee~ £/ ~ e 3 , leading to a.\ = 3. 

Consider now a type-3 bifurcation that occurs at some 
A 3 . We write A = A 3 + e and adapt our previous argu- 
ment. Let m(t) be the mass of the newly-formed central 
cluster, and let M be the final mass of the two major 
clusters surrounding it. We have m = —mM, since the 
central cluster interacts with half of the mass M on either 
sides. Therefore 



m(t) 



-Mt 



(9) 



In contrast to Eq. (g) where the amplitude was of or- 
der e, the amplitude in Eq. (^|) is of order unity. This 
arises because the range of opinions that contributes to 
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the ultimate central cluster is of the order of the interac- 
tion range. Now the argument proceeds as before. The 
width of the large cluster varies as e~ Mt / 4 . The condi- 
tion for the central cluster and its neighbors to decouple 
is e ~ Mt // 4 ^ e . At this point, we have m/ ~ e~~ Mtf ~ e 4 , 
resulting in the exponent 0:3 = 4. 

The heuristic argument we have presented is consis- 
tent with the extremely small mass of the minor clusters. 
For large A, the system is governed by the parameter 
1=\L — 1 « 0.08, the excess between the adjacent clus- 
ter separation and the interaction range. This parameter 
essentially plays the role of e, the small distance from a 
bifurcation point. The two extreme minor clusters evolve 
according to the mechanism that led to the e 3 behavior 
near a type-1 bifurcation. Thus, their mass can be es- 
timated m! ~ e 3 w 5 • 10~ 4 . On the other hand, minor 
clusters in the bulk evolve according to the mechanism 
that led to the e 4 behavior near a type-3 bifurcation. Ac- 
cordingly, their mass is estimated by m ~ e 4 w 4 • 10~ 5 , 
again a reasonable value. 



isfy two conservation laws: 

N N 

J2 p i = 1 > E iP * = A ' ( n ) 

»=1 <=1 

with 1 < A < N. The former (latter) reflects conserva- 
tion of the total population (opinion). As a result, there 
are A — 2 independent variables for an A-state system. 



A. Typical Behavior 

Equations (|l0|) are non-linear and therefore for A > 4 
they cannot be solved to obtain explicit formulae for 
Pi(t). However, the qualitative behavior can be still un- 
derstood. For example, Eqs. ( |Io|) admit only the simplest 
type of attractors - fixed points - while limit cycles are 
impossible. We illustrate this by analyzing small values 
of A to highlight the new qualitative features that arise 
as A increases. 



III. THE DISCRETE VERSION 

Often, one faces a choice among a finite set of options, 
so it is natural to consider a discrete version of the com- 
promise model. While interesting on its own, the discrete 
model also enables us to illuminate many qualitative as- 
pects of the behavior in the continuum case. Discrete 
systems are governed by a finite set of non-linear rate 
equations, so explicit solutions are generally impossible. 
Nevertheless, we can gain considerable insight by investi- 
gating small systems, using stability analysis and related 
tools from theory of ordinary differential equations p7| . 

In the discrete version, each agent can take on an opin- 
ion from a set of A equally-spaced values. To impose an 
interaction threshold and also to ensure that the out- 
come of an interaction remains within the state space, 
two agents interact as follows: (a) If the opinion differ- 
ence is greater than two, there is no interaction, (b) If 
the difference equals two, the agents reach a fair compro- 
mise and each takes on the average opinion value, (c) If 
the opinion difference equals one, nothing happens. 

We label the opinion states as i — 1,2,..., A, so 
schematically, in a compromise event (i— 1, i+ 1) — * (*,«)■ 
Denote by Pi (t) the fraction of the population that has 
opinion state i at time t. For general N the fractions 
Pi(t) obey the rate equations 

P = 2P,_iP i+ i - P(Pi- 2 + P+2). (10) 

This equation formally applies for i at least 2 spac- 
ings away from the boundaries (at 1 and A). Setting 
P_i = P = Pv+i = Pv +2 = in Eq. Q, yields the 
governing equations near the boundaries: Pi = — P1P3, 
P 2 = 2PP 3 - P 2 P 4 , Pat = -PvPv-2, and Pv-i = 
2P/vP/v_2 — P/v-iP/v-3- Again, the fractions Pi{t) sat- 



1. Isolated fixed points 

For A = 3, there is a single fixed point located at 

f (2- A, A- 1,0) whenl<A<2, 
\ (0,3 -A, A -2) when 2 < ,4 < 3. 

This point is stable. Asymptotically, it is approached 
exponentially fast in time; e.g., for 1 < A < 2 one finds 
P3 oc e^ 2 "' 4 )'. An exception arises for the symmetric 
initial condition (A = 2) when the final central state 
(0, 1,0) is approached algebraically in time: Pi = P3 — > 
t~\ 

For N = 4, there is also a single stable fixed point 
located at 

(2- A, A- 1,0,0) whenl<yl<2, 
(0,3- A, A- 2, 0) when 2 < A < 3, 
(0,0,4- A, A -3) when 3 < A < 4. 

Additionally, there is a fixed point ( 4 -^ 4 -,0, 0, ^jp-) that 
is always unstable. The stable fixed point is approached 
exponentially in time. 



2. Lines of fixed points 

For N = 5, some fixed points are no longer isolated 
but instead they form lines. Indeed, Eqs. ( |Io| ) admit 
fixed points of the generic forms (P* , P 2 *, 0, 0, P 5 *) and 
(P*, 0, 0, P 4 *, P 5 *) that are stable when, respectively, P.* > 
3P 5 * or P 5 * > 3P*. Recalling the conservation laws ( |Tl| ) 
we can write these fixed points in the form 

/ (2 - A + 3P 5 *, A - 1 - 4P 5 *, 0, 0, P 5 *), , , 
\ (Pf*, 0, 0, 5 - A - 4P^, A - 4 + 3P^) y ' 
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and obviously the fixed points form lines 
points 



The fixed 



2|) from the first line are stable when 1 + 4P 5 * < 



A < 2, the fixed points (12) from the second line are sta- 
ble when 4 < A < 5 — 4P* . There is also a stable isolated 
fixed point located at 



(0,3- A, A- 2,0,0) 
(0,0, 4 -A, A -3,0) 



when 2 < A < 3, 
when 3 < A < 4. 



(13) 



For 2 < A < 4, every initial condition is in the basin 
of attraction of the isolated fixed point ([l3"|). There- 
fore we know the fate of the system without explicitly 
solving the rate equations. This statement tacitly as- 
sumes that a trajectory does not approach a limit cy- 
cle or other complicated attractor; this will be justified 
later. In the complementary range A G (1,2) U (3,4) 
the trajectories approach one of the stable fixed points 
in (|l2]). For example, if 1 < A < 2, the final state 
is (2 - A + 3P 5 *,A - 1 - 4P*,0,0,P 5 *) for some P* G 
(0,^7^)- To determine which fixed point is actually 
reached depends not only on the initial average opinion 
A = iPi (0) but also on other details of the initial con- 
dition and requires a complete solution. Qualitatively, for 
an initial state that is central in character (2 < A < 4), 
the final occupation fractions are concentrated in a sin- 
gle central cluster consisting of two adjacent sites. Con- 
versely, for an initial state that is biased toward one ex- 
treme, the final state consists of two extremal clusters. 

The above elementary examples demonstrate a simple 
principle. The rate equations have multiple stable fixed 
points. Each stable fixed point is a basin of attraction for 
some region in the space of initial conditions. The dy- 
namics determine which stable fixed point is eventually 
approached. In the continuous version, a similar situ- 
ation occurs where there are enormously many steady 
states of the form (||). Moreover, we see how depending 
on the initial conditions, the system can reach a single 
central cluster or two off-center clusters. 

In the remainder of this subsection we consider sym- 
metric situations, Pi = P/v+i-i- For an TV-state system, 
we can choose 1,2,..., \N/2~\ independent states, where 
fJV/2"| is the smallest integer that is larger than or equal 
to N/2. Conservation of population diminishes the num- 
ber of independent variables by one, while the second 
conservation law is redundant as A — (N + l)/2 always. 



3. Explicitly solvable case 



For N = 6, the rate equations (|10|) are exactly solvable 
and the solution neatly illustrates the features described 
in the previous subsection. Using symmetry and normal- 
ization, we can treat the system in the two-dimensional 
triangular domain defined by (Fig. |^) 

T = {(Pi,P 2 )|Pi > 0,P 2 > 0,Pi +P 2 < 1}. (14) 




FIG. 6: Schematic (Pi,P2) phase plane for the symmetric 6- 
opinion system. Shown are the isolated fixed point (dot), the 
line of fixed points (heavy line - dashed for unstable and solid 
for stable), and the separatrix (dotted) that demarcates the 
basins of attractions of these two sets. 



P 2 * = 1. Linearizing near the isolated fixed point we find 



that P = (Pi, Pa)"* satisfies 



with M 



-1 
2 -1 



(15) 



from which the origin is a degenerate stable node. Lin- 
earizing near a fixed point (P*,P 2 *) we find that it is 
stable iff Pf > P 2 * and unstable iff Pf < P 2 *. Thus the 
isolated fixed point has a finite-area basin of attraction, 
while every point (PfjPf) with Pj* > P 2 * has a basin of 
attraction that is a one-dimensional manifold (Fig. |^). 
In principle, a two-dimensional system could have closed 
orbits. However, every closed orbit in a two-dimensional 
system must enclose fixed points p7|| . Here, all fixed 
points lie on the boundary of the phase plane T, so closed 
orbits that encircle a fixed point are impossible, thereby 
ruling out cycles. 

The solution to Eqs. ( [T~0| ) for N = 6 with symmetric 
initial conditions is given in the Appendix. This solution 
gives the following behavior i n th e phase space T (see 
Fig. 5). There is a separatrix (A6) that joins (P 1; P 2 ) = 
(1/2, 1/2) with (\/e/4, 0). The part of the phase plane to 
the left of the separatrix is the basin of attraction of the 
isolated fixed point, while the complementary region is 
the basin of attraction of the line of fixed points (Pj*, P 2 ) 
with P* > | > P 2 . These fixed points are approached 
exponentially in time. Finally, the separatrix itself is the 
basin of attraction of the fixed point (1/2, 1/2). In this 
borderline case the relaxation is algebraic rather than 
exponential: 



Pi 



2t~ 



(16) 



There are two kinds of fixed points: An isolated fixed 
point (0, 0) and a line of fixed points (P*, P 2 ) with P* + 



Both consensus and polarization are possible outcomes - 
which actually occurs depends on the initial condition. 
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4- Large N 

For N > 7, the systems are \N/2] — 1 > 3 dimen- 
sional, and already the trajectories of 3-dimensional sys- 
tems may exhibit a vast range of behaviors including 
chaos |27], [2S|| . In the present case, however, we find 
that there are simply more and more fixed points, and 
they appear as isolated fixed points, lines, surfaces, and 
higher-dimensional sub-manifolds. 

For N — 7, the system is three-dimensional, and the 
phase space is the simplex 

S = {(Pi,P 2 ,P 3 )\P j >0,Pi + P 2 + P 3 < 1}. (17) 

For simplicity, we denote states by (P 1; P 2 , P3, P4). The 
system admits the following fixed points: 

1. A line of fixed points (P-j^P^, 0,0) corresponding 
to a polarized society. 

2. A line of fixed points (P^OjO, P4) corresponding 
to a society with both centrists and extremists. 

The second case includes the central consensus state 
(0, 0, 0, 1) as a special case. 

Linearizing around the fixed point (P^P^, 0,0) we 
find that P = (P3,Pi) T satisfies 



M 



-PI 2P 2 * 
-2P,* 



(18) 



implying that (P^Pj, 0,0) is a stable node (we tac- 
itly assume that P X *,P 2 * > 0; when Pj* = 2P 2 * this 
node is degenerate). Linearizing around the fixed point 
(Pi*, 0, 0, P|) we find that P = (P 2 , P 3 ) T satisfies 



dP 

dt 



MP, 



M = 



-P 4 * 2P{ 
2P 4 * -PT 



(19) 



implying that (P*,0,0, P4) is a saddle point. Therefore 
the fixed points (P*, 0, 0, P4 ) are unstable (again it is as- 
sumed that P*,P% > 0). The two extreme fixed points 
(1/2, 0, 0, 0) and (0, 0, 0, 1) are neutrally stable in the lin- 
ear approximation. Therefore one must go beyond the 
linear approximation to probe the stability of consensus. 
Numerically, one typically finds that the system reaches 
consensus {e.g., starting from the uniform initial condi- 
tion). Therefore, consensus appears to be stable. 

For larger N, we determined the final state numeri- 
cally. To compare with the continuum case, we start with 
the uniform initial condition. Generally, the final state 
consists of non-interacting clusters. Each cluster consists 
of a pair of occupied sites and clusters are separated by 
at least two empty sites. We assign each cluster a mass 
m equal to the combined occupation of the two sites, and 
a position x determined from a weighted average. 

As a function of N , the number of clusters grows via 
a series of transitions, rather than bifurcations (Fig. [?]). 
The main difference with the continuum case is that while 
minor clusters occasionally appear, they do not persist 
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FIG. 7: The location of the final clusters in the discrete model 
as a function of TV, with TV odd. 



in the form of minor branches. Otherwise, there are 
many similarities. Transitions involving major and cen- 
tral clusters are observed; in particular, there are type- 
2 ({0} — > {— x,x}) and type-3 (0 — > {0}) transitions. 
These transitions are arranged in a periodic structure 
2,3,2,3, and the transition diagram is approximately 
invariant under the transformation N — » N + Nq with 
TVo = 12. Branches of major clusters carry almost equal 
masses, and remarkably, despite the discreteness, these 
branches grow linearly with N. 



B. General Features 

1. Volume Contraction 

The system of rate equations is dissipative, that is, vol- 
umes in phase space contract under the flow. Generally 
for a system of differential equations Pj 
V(t) changes according to 



Fj, a volume 



dV 
~dt 



F-ndS 



V -FdV 



(20) 



where n is the outward normal on the bounding surface 
S(t) that encloses V(t) and V • F = For thc 

OF- 

infinite discrete system -^p- = —Pj-2 — Pj+2 and thus 
the contraction rate is twice the (conserved) total pop- 
ulation: — Y] . -gp- — 2j^ fe Pfe. Hence volumes in phase 
space shrink exponentially in time. If we set the total 
population equal to 1, we get V(t) = V(0) e~ 2t . 

For finite systems, the contraction rate is generally not 
a constant but nevertheless volume still shrinks exponen- 
tially in time according to the bounds 



V(0) e~ M < V(t) < 1/(0) e 



(21) 
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For example, for N = 4 the contraction rate is J^k 
I : therefore, V(t) — V(0) e~*. For N = 5, the contraction 
rate is I+P3; for N = 6, the contraction rate is I+P3+P4, 
ete. This is consistent with the system evolving toward 
fixed points. Nevertheless, cycles and strange attractors 
are possible in volume contracting systems, with the cel- 
ebrated Lorenz system being a prime example . 

2. Lyapunov functions 

We now demonstrate that our system has only fixed 
points (many of which are actually fixed submanifolds 
in phase space) by constructing a Lyapunov function 
L = L[P(t)], viz. a smooth function that decreases along 
trajectories. The existence of a Lyapunov function rules 
out cycles. Indeed, suppose that there is a periodic so- 
lution with period T, then the integral dt ^ over the 
period would be negative because the Lyapunov function 
is decreasing. On the other hand, the integral must be 
equal to zero since the trajectory returns to the starting 
point. This contradiction means that no periodic solu- 
tions can exist. 

Consider, e.g., symmetric situations. For N = 7, 

L = P 4 + 2P 3 + 4P 2 + 8P1 (22) 
is a Lyapunov function; indeed, it satisfies 

^ = -2P 1 P 3 -2P 2 P 4 , (23) 
at 

so the derivative is strictly negative inside the simplex S 
(it vanishes only on the two lines of fixed points on the 
boundary of S). 

Generally, we can construct Lyapunov functions for all 
N. For instance, when N is odd we write N — 2M — 1 
and verify that 

M 

L = ^ M - j Pj, (24) 
3=1 

is a Lyapunov function as it satisfies 

, M-3 

d ± = -2P M _ 2 P M Yl 2" ' ; 'PP.,- (25) 



3. Negative Diffusion Instability 

In the absence of boundaries, any uniform state is a 
trivial solution of the nonlinear set of the ordinary differ- 
ential equations ( |io| ) . To check the stability of the uni- 
form state, Pi = const., we treat i = x as a continuum 
variable. Writing P(x,t) = 1 + <j>{x,t) with (f)(x,t) <C 1, 
this perturbation evolves according to 

<t>t+[<t>+l<t>x X + l<t> 2 ^ -0, (26) 



where the subscripts denote partial differentiations. To 
lowest order, this is the diffusion equation with a nega- 
tive diffusion coefficient. Hence, the uniform state is un- 
stable to the perturbation <fr(x,t) — exp(ikx + Xt) when 
k < y/6/7. Therefore, minute details of the initial condi- 
tions are magnified, ultimately resulting in isolated clus- 
ters. However, the nonlinear terms in Eq. (|2^) eventually 
counter the instability. 



IV. DISCUSSION 

The interplay between compromise and conviction 
leads to intriguing opinion dynamics. The system ul- 
timately reaches a static state that consists of a finite 
number of noninteracting opinion clusters, and the num- 
ber of these clusters increases via an infinite sequence of 
self-similar bifurcations as the opinion range increases. 
In the bulk of the system, clusters are organized in a pe- 
riodic lattice of alternating minor and major clusters. A 
central cluster may or may not exist, and its size exhibits 
a complex periodic behavior. 

As a model of mathematical sociology, the compromise 
model is appealing in its simplicity, yet its behavior is fa- 
miliar in everyday experience. A political system may 
or may not contain a centrist party. Alternatively, it 
may consist of two (or more) off-center parties. Further- 
more, the existence of marginal parties halfway between 
two major ones is also reasonable. Artificial features of 
the model, such as the identical separation between par- 
ties, can be easily circumvented by introducing hetero- 
geneities. For example, since different agents may have 
different levels of conviction, it may be natural to have 
interaction thresholds that are specific to each individual. 

As a dynamical system, the compromise model exhibits 
the simplest types of attractors, namely, fixed points that 
are either sinks or saddles. In the discrete case, we con- 
structed Lyapunov functions and also established that 
limit cycles and strange attractors are impossible. This 
conclusion extends to the continuum case. The second 
moment decreases monotonically with time and hence, 
it is a Lyapunov functional. Generally, each stable fixed 
point is the basin of attraction of some region in the space 
of initial conditions. In other words, the rate equations 
map an initial state into a final state. Given the large 
number of these states, it is not obvious how to charac- 
terize such a map. One practical approach is to obtain 
statistical properties of the final state by averaging over 
all possible initial conditions. In the discrete case, we 
find that starting from a random initial state ({P;(0)} 
randomly chosen in the A-dimensional hypercube) the 
distribution of cluster masses and the distribution of the 
separations between clusters in the final state are inde- 
pendent of the system size, for large enough systems. 

There are important questions concerning robustness 
of the bifurcation diagram with respect to variations in 
the dynamical rules or in the initial conditions. We ex- 
amined only the former. When the opinion difference 
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between two agents is merely reduced by a fixed factor, 
i.e., they reach partial compromise, an almost identical 
bifurcation diagram is found. The effects of asymmetry 
in the initial conditions, and the dependence of the lo- 
cation of bifurcation points on the shape of the initial 
distribution deserve a careful study. 

There are numerous possible generalizations of the 
compromise model. One is to increase the dimension 
of the opinion space. Do the final opinions form a lat- 
tice? and if yes of what type? Yet another open question 
is the role of spatial dimension. In the present work, 
we implemented the mean-field limit where any agents 
equally well interact with any other agent. However, if 
agents are located at lattice sites with interactions only 
between nearest neighbors, spatial correlations are ex- 
pected to emerge |2Jj, |30| . We anticipate that for the 
discrete system in which N > N c (d), with N c (d) depend- 
ing on the spatial dimension d, the system will freeze into 
a large number of non-interacting opinion domains. The 
case of continuum opinions is unexplored, but we expect 
both slow dynamics and coarsening of opinion patterns 
as the final state is approached. 
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Depending on the initial conditions, P;(t) either re- 
mains positive or it reaches zero. In the former case, 
P?i{ T ) — ^ 1 as t — ► oo and asymptotically the isolated 
fixed point is reached. To determine the approach to the 
fixed point in terms of original time variable we write 



t = 



dr' 



lo W 

with P3 given by (|A2|). Asymptotically we find 

t = t-c+0{te~*), 

[2P 1 (0)t + P 1 (0) + P 2 (0)] e~ T 



(A3) 



dr 



1-[2P 1 (0)t + P 1 (0) + P 2 (0)} 



Substituting this into (A.2) we arrive at 



-2t\ 



Pi{t) = nie-* + C(te 
P 2 (t) = [2IIi t + n 2 ] e~* + G(t 2 e 



-2t\ 



with III = Pi(0) e c , and Tl 2 = [P 2 (0) - 2cP 1 (0)] e c . 

In the complementary situation, P^(t*) — at some 
t* and then it becomes negative. I n th e limit r — > r* 
the physical time diverges; see Eq. ( |A3| ). Therefore the 
range r > t* is physically forbidden so that the system 
reaches a fixed point (P*, P 2 *), with 

P; = Pi(0)e- T \ 

P; = [2P 1 (0)r*+P 2 (0)] e- T ". 



APPENDIX A: SOLUTION TO THE RATE 
EQUATIONS FOR N = 6 

The symmetric N — 6 system simplifies after replace- 
ment of the time variable t by 



T= [ dt'P 3 (t'), 
JO 



(Al) 



as the rate equations reduce to the linear system 
P[ = -Pu P 2 = 2P 1 -P 2 , Pi = P 2 -Pi, 

where ' = djdr. Solving these equations we obtain 
Px(r) = P!(0)e- r , 

P 2 (r) = [2P 1 (0)T + P 2 (0)]e- T , (A2) 
P 3 (r) = l-[2i\(0)T + i\(0) + flj(0)] e- T . 



The approach to this fixed point is exponential in time. 

The borderline between these two regimes occurs when 
p3( T ) > for all t r* , i.e., the curve Ps(r) touches the 
r axis horizontally. Thus we require both 



and 



I = [2P X (0) r* + P 1 (0) + P 2 (0)] e~ T * ( A4) 



2Pi(0) = 2P 1 (0)r*+P 1 (0) + P 2 (0). (A5) 



The second relation gives t* = [Pi(0) - P 2 (0)]/[2Pi(0)]. 
Substituting this into (A.4) yields the separatrix 



Pi(0)-P 2 (0) 
2Pi(0) 



= ln2Pi(0). 



(A6) 
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